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The well known string solution to Bethe ansatz equations is shown to be inconsistent in its 
widely accepted form. A valid refinement demanding higher order corrections in subsequent roots 
is identified. A new alternative string solution is proposed for finite long chains, consistent with the 
non-Bethe strings observed earlier. 

PACS numbers: 03.65. Ge , 75.10.Jm , 05.90+m , 11.55.Ds 
Seventy years ego a phenomenal paper by Hans Bethe jj], opened up a new direction in physics, namely the theory 
of exactly solvable quantum systems. The pioneering concepts of coordinate Bethe ansatz (CBA) and string solution 
■ to Bethe ansatz equations (BAE) at the thermodynamic limit, introduced in this paper for the isotropic spin-i chain, 
proved to be an universal tool for almost all such models discovered in subsequent years. Next development, though 
^ • came rather late, solved a number of important problems using the same method. Some of such milestone works, 
where a wide range of one dimensional models were exactly solved at the finite chain as well as in the thermodynamic 
limit may be listed as follows. Bose gas interacting through 5-function potential ||]; XX Z-sp'm chain || along with 
the application of thermodynamic Bethe ansatz (TBA) Q; repulsive 5- function fermion gas MM; ID Hubbard model 
£N) ' @; fully anisotropic XYZ spin chain Q etc. In most of these models the Bethe string plays a crucial role, especially 
in the TBA treatment. A systematic application of TBA to a variety of models has been considered in details in ||. 
Another remarkable achievement is the algebraic formulation of Bethe ansatz (ABA) 10 Ol w hich is more suitable 
for quantum field models like nonlinear Schrodinger equation (NLS) |0|, sine-Gordon (SG) |13[| and Liouville models 
|l4j |. Though ABA differs considerably from its coordinate formulation (CBA), the end equations, e.g. BAE turn 
\ out to be equivalent. Therefore the Bethe string solutions may appear also in such integrable field models forming 
quantum solitons or breathers ]l0| , ^3[ . The widely accepted and frequently quoted form of the Bethe string is given 
g i by HQ 

\i=\ + i(l- 1 -^)+iO(e- aN ), l = l,...,r, (1) 
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with a positive, which is supposed to satisfy the Bethe ansatz equations 

/ \ N 

hill) = TT ^-A w +i I=1 (2) 



for large N. It is evident from the above solution that exponentially small corrections: 0(e~ aN ) are needed for large 
but finite values of N. Moreover, the form (|l|) suggests that these corrections are of the same order for all the roots, 
i.e. they should be at least as Sie~ aN . Consequently, in the first approximation, i.e. for N = oo, all corrections must 



C$ ' vanish making the string solution exact and given in the obvious form |lG] |l8| 

A, =\ + i(l-Z±l), 1=1,. ..,r (3) 

We however find surprisingly that both the above well known and well accepted forms (Q) and ([}]) are, strictly 
speaking, not consistent with the BAE (||) for higher r (r > 2) values. We observe further that a fine-tuned form with 
carefully defined approximations, where higher and higher order corrections appear in the subsequent roots, can only 
pass for a valid solution. Such inconsistency in the commonly accepted form of the string solution and the necessity 
for its refinement perhaps have never been addressed before. On the other hand, in early eighties a series of papers 
reported results where for long but finite chains some non-Bethe string solutions were found to appear. 

We intend to propose here a new string solution also for finite long chains, which is consistent with the BAE without 
finely adjusted orders of correction. Roots within such strings may group together in sets of two, three and four 
in the thermodynamic limit, providing an explanation for the similar non-Bethe string structures observed earlier 
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|_5|,^6|Jl9| . An important feature of this alternative string is that, the spacings between the imaginary parts of its 
roots, in contrast to the Bethe string, are no longer same for all solutions but depends on the spin excitation number 
r. As a consequence the string length does not grow linearly with r as in the standard case, but saturates as 
giving only short strings. 

For concreteness we consider the Heisenberg ferromagnet with nearest neighbor interactions: H s = — ^ 2^2n=i ' 
(?„ + i . We look briefly into the Bethe ansatz method for solving its eigenvalue problem to see how BAE (|^) 
arise, before concentrating on the string solution (|l|). The Bethe ansatz for the eigenfunction in r-down spin state: 
^ r = J2{nn} a ( m i! m 2i ■ ■ ■ j m r )\mi,m 2 , ■ ■ ■ , m r > is chosen as 

r 1 

a(m 1 ,m 2 , . . .,m r ) = ^%xpz[^ fc P; ra ; + - <M,pJ, (4) 

P 1 = 1 Kn 

where P is any permutation of the numbers 1, 2, . . . , r and Pi is the number that replaces I under this permutation. 
It is easily seen that if none of these spin excitations are adjacent, they become noninteracting and the Schrodinger 
equation H^ r = E r ^> r leads to 

r 

E r a(mi,m 2 , ■ ■ ■ , m r ) = . . , rrij + 1, . . .) + a(. . . , rrij — 1, . . .) — 2a(. ..,m,j,...). (5) 

i=i 

The ansatz (|j) representing a superposition of plane waves is clearly an exact solution of (^) for arbitrary phases. 
However as soon as any two of the spins become adjacent, say m; + i = mi + 1, interaction sets in making the 
equation different from f^). These two sets of equations can only become consistent if the condition on the amplitude 
2a(. ..,mi, mi+i, . . .) = a(. . . , mj, mi, . . .) + a(. . . , mi + 1, mi + 1, . . .). is satisfied. This in turn fixes the phase 
through pseudomomenta as cot<f)i n = ^(cot^ — cot^f), —n< <f>i n < +7r with the obvious antisymmetry property 
4> n i — —4>in- The next problem is to identify the equations for determining the variables ki, I = 1, 2, . . . , r, for which 
we use the periodic boundary condition a{m\, m 2 , ■ • • , m r ) = a(m 2 , ■ • ■ , Tn r , mi + N) and involve ansatz (^) to get 
an important set of equations 

e *N = e 'E„^- l = l,...,r, (6) 
known as BAE. Introducing rapidity variables as 

\ 1 M • • i<h ln Xi-Xn+i 

x i = o cot T> ^ving e = t ; : (7) 

2 2 A; — A„ - I 

one can easily rewrite BAE (||) in the form (^). It is worth noting that in the ABA method the BAE appear 
automatically in the form (||), where A;'s correspond to fixed values of the spectral parameter. The solutions for Az 
to BAE (Q) must be equivalent to the solutions fa for (Q), though their ranges are different. The total number of 
solutions of BAE, as shown in [jjj, must be given by real as well as complex roots. To find the complex solutions for 
fa = ui + ivi , we notice that for positive (negative) values of vi the LHS of (||) exponentially vanishes (diverges) for 
large N. Therefore at least one of the 4>i n 's must have positive (negative) imaginary part dominatingly large as O(N). 
The well known 7'-Bethe string solution, which is supposed to fulfill this criterion is given in the form ([!]) or (||) for 
variables A;, as mentioned above. 

Let us however check the validity of this solution for all values of r by direct insertion into BAE (||) . In the simplest 
case of r = 2 we get from (JsJ) in the first approximation: Xi = X — i^, X 2 = A*, forming a complex conjugate pair. 

Therefore it is enough to check Eqn (^) only for / = 1, where LHS = f ^-i ) gi yes a vanishing term at N = oo. This 

is easily seen by taking its modulus. RHS = on the other hand vanishes due to Ai — A2 = —i, proving thus 

the 2-string as a valid solution. For r = 3 with A2 = Ao real and Ai = A3 = Ao — i, we have to check equations for 
1 = 1,2. Here again for the neighboring pairs Ai — X 2 = X 2 — A3 = —i and consequently, the equation for I = 1 holds 
on similar grounds as above at large N. For I = 2 on the other hand, zeros appearing in the RHS from the terms with 
n = 1 and n = 3 get canceled giving a nonvanishing term. This however does not lead to any contradiction, since A2 
being real the LHS is a pure phase and does not vanish for any N. 

Since the Bethe string in the form (Q) is valid, as shown above, for r = 2 and r = 3, it is probably natural to 
assume that it will be valid for arbitrary r, by analogy. This assumption is generally made in most of the works in 
proving Bethe string solution for higher r p^JT^ ] . We show however that this assumption is not true in general and 
for demonstration consider the case r = 4 in some detail, where the Bethe string obtained from (^) is of the form 
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Ai = X a -i^+iO(e- aN ) : A 2 = A - i~ + iO( e - aN ), with A 4 = \{, A 3 - A*, (8) 

and one has to check only for the independent solutions with I = 1, 2. Defining for convenience a function V s (x) — 
with the obvious property 

Vl ( T i + 0{e- aN )) = 0{e^ aN ), (9) 

we see for I = 1 that the LHS of (|): V^Ai) = (x^ffH vanishes exponentially for large N, while the RHS 

is a product of three factors Vi(Xi — A 2 )Vi(Ai — As)Vi(Ai — A 4 ). Since only neighboring rapidities can give Ai — 
A2 = — i + 0(e~ aN ), we conclude immediately that the RHS also vanishes exponentially due to the single factor 
Vi(Ai — A2) = 0(e~ aN ) and the remaining giving finite contributions. This clears the string solution for I = 1. 

However, as we see below, the contradiction starts from 1 = 2,. The LHS in this case: Vj^^) = f A ^ j is again 

exponentially small and the RHS contains three factors Vi(A2— Ai)Vi(A2 — \3)Vi(\2 — A4). However the important point 
is that, there are now two neighboring roots of A2 contributing as A2 — A3 = — i + 0(e~ aN ) and A2 — Ai =i + 0(e~ aN ) 
and consequently using property (g) we get Vi(A2 — A3) = 0(e~ aN ), while Vi(A 2 — Ai) = 0(e aN ) and Vi(A 2 — A4) as a 
finite term. Therefore the important conclusion we arrive at is that, all singularities in the terms of RHS get mutually 
canceled leaving only a finite contribution even at N — ► 00, contrasting the LHS. It is apparent now from the above 
reasoning that this inconsistency would always appear for higher r, in the Bethe string solution in its accepted form 
(|l|) for all its roots A;, with I = 2, . . . , r — 1, having two nearest neighbors, giving 

LHS = VFiXt) ~ 0{e~ VlN ), (10) 

2 

while RHS~^(A,-A, +1 )Vi(Ai-A,_i) = 0( e - aN )0(e aN ) ~ 0(1). (11) 

Only for the end roots having single neighbors or for roots with real values, the form (|]J) remains consistent. Since 
for r — 2, 3 all roots fall into this category, they give valid string solutions. 

The next natural question is how to refine the structure of (B), so that we can get valid Bethe string for all values 
of spin excitations r. For answering this question we see from (ll]) that the reason for failure of the string form (|l|) 
is that, it allows the same order of correction 0(e~ aN ) for all its roots A;. Therefore, the refined form for the Bethe 
string, consistent for arbitrary r should be given as 

\ l= \ Q + i (l~?-t±) +l O(e- a ' N ), and A r+ i_; = A,*, / = l,...,s, (12) 

with r = 2s for even and r = 2s + 1 for odd r with real X s +i = Xq. The essential point for the validity of this refined 
string is that, the exponential orders in ( [l2| ) must be fine-tuned as a strictly growing sequence 

< ai < • ■ ■ < ai < ai + i < • • • < a s , (13) 

and more precisely they should be given by the recurrence relation 

ai-a l ^ 1 =v h where v t = - Zn(l + Ki) > 0, (14) 

with K[ = A ^ 1 r ~_^^ 2 > for I = 1, . . . , s. To check the validity of the refined form ([l^) we insert it again in BE A (^) 



and notice that the LHS ~ 0(e Vl ) at large N similar to (|10| ), while the RHS is now given by 



RHS ~ V 1 (X l - A; + i)Vi(A ; - A/_i) = {0( e -^ N ) - 0( e - a ^ N )) (0( e -^ N ) - 0( e - a ^ N )) 1 

w 0{e- aiN )0{e a '- lN ) - 0{e- VlN ), (15) 



{e~ a ' N ) 



which proves the claim. In deriving the expressions in ( [15] ) we have used the relation ([wj) and approximated 0(e 
0( e -°"+ lN w 0(e~°" N ) and 0{e~°" N ) - 0{e~°"- lN w 0(e~ a '- lN ) neglecting higher order smaller terms due to strict 
inequalities (|T3|). It is important to note that since all higher orders of correction are present simultaneously in (|12j), 
it is not principally possible to approximate this solution order by order for r > 4 cases. For example, one may 
naively consider the first approximation by putting all corrections 0(e~ a ' N ) = in (12) reducing it to the form 



(|3J), which is however not a valid solution for r > 4, as we have seen above. In any approximation, one can at best 
neglect the highest order of correction by putting 0(e~ QsN ) = 0, keeping however all other corrections as in (pj). 
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The limit N = oo should be inserted in the solution only at the end after performing all calculations, even in the first 
approximation. For r = 2, 3 however, since 8 = 1, we can neglect all corrections in ( |l2| ) in the first approximation and 
get (||) as the consistent solution. In fact for these cases only the form ( |T^) coincides with (|l|). 

It is fair to mention here that, though the Bethe string in the form ( |l2j ) refined through jl^), ([Ti]), as far as we 
know, has not appeared in any of the earlier works, the form presented by Gaudin Q agrees qualitatively with (p^), 
though the requirement (jl^), ( |l4| ) necessary for the validity of the solution was not specified in his paper. Note that 
the complex string solution for pseudomomentum k[ equivalent to (Jl^) should be given through the mapping ([?]), 
from which it is clear that in the corresponding phase iJ2n^i 0i« appearing in (|^), the terms 4>u±i must have large 
imaginary parts having corrections like iftiN and —i(3i-\N, with the inequalities 

0< A <•••<$_!< A <■■■<&, Pi -0i-i= vi, (16) 

reflecting @, (0)- Bethe [| suggested also that at least one of the 4>in has to have a very large imaginary part of 
O(N). However the essential requirement ([l6|) for the consistent higher r strings was not emphasized and @) was 
taken as the solution for arbitrary r in jjj (see Eqn (30)). We have shown however that one can not get any solution in 
the form (|§|) for r > 4, since the corrections can not be neglected from the beginning, even in the first approximation. 

Thus we have identified a consistent but complicated form ([l2| ) for the Bethe string, which demands simultaneous 
involvement of all higher order corrections (up to s) for its validity. However such fine-tuned structures growing 
linearly with r might be on physical grounds less probable to appear in low lying excitations. Earlier findings of short 
strings with non-Bcthe nature seem to support such arguments. This therefore motivates us to find some new type 
of string solution with desirable properties like shorter string length; possibility of order by order approximation and 
finally possible explanation of non-Bethe strings observed earlier. The alternative string we propose has the form 

A, =\v+i-(l- 1 ^-)+iO(e- aN ), l = l,...,r. (17) 
s 2 

Note that the corrections appearing in all its roots have the same order and the spacing between the neighboring roots 
d(l) = j varies inversely with r. In contrast to the Bethe string with constant spacing d(l) = i, the string jl7| ) with 
more roots is more dense along the imaginary axis and its length: L{r) = 2( -' * ^ , measuring the distance between end 
roots, is bounded between 1 > L(r) > 2. Recall that for Bethe string L(r) = r — 1. 

For checking the validity of ( JlTf ) we put it directly in BAE (||) and observe that while the LHS is similar to (|Io| ) 
giving exponentially vanishing term Vt'(Xi) for large N, the RHS behaves differently from the standard case (|ll|). 

2 _ _ 

Since now d(s) = i, the pairing partner having the crucial difference A; — A; = — i + 0(e~ aN ) is given by A; = \i+ s - It is 
also seen easily that in this case for every complex root A; such pairing can occur only once, i.e. there can be no other 
root except A; giving either i or —%. This is in sharp contrast with the standard Bethe string (Q), where singularities 
can come from two neighboring roots. Therefore unlike jll|), we get for ( |l7| ) RHS ~ Vx(Xi — A; +s ) = 0(e~ aN ), i.e. 
an exponentially vanishing term like LHS. It is also worth noting that since the alternative string ( |l7| ) allows order by 
order approximation one can consider the first approximation by putting N — oo, when the correction term 0{e~ aN ) 
drops out from all roots of (|l7j), making it an exact solution. Recall again that this is not possible for a consistent 
Bethe string given in the form (|l2|). 

From the above arguments it is clear that for each A; in ( p"^ ) its partner A; = \i+ s and its conjugate — A r +i_; 
are also included in the solution and since the partner of its partner and conjugate of its conjugate is the root itself, 
the group of four {A;, A;, A*, A ; *} forms a close unit, being the decisive contributors to the equation at large N and 
thus becomes almost independent entries. Therefore all the r roots can be thought of to be broken up into units of 
four at the thermodynamic limit, providing an explanation of the quartet formations found earlier. However when 
the number r is not divisible by 4, doublet and triplet groups can also occur. For I = ^-j^ when r even and for 
I = when r odd, one gets A; = A ; *, and consequently it can group only as a unit of two, forming a two string. 
For odd r, A s +i = Ai = is real and gives a triplet {Ai,Ai,A*}, forming a 3-string. So the general rule is that, 
in the thermodynamic limit r = An results n quartets, r = An + 1 gives n — 1 quartets, a doublet and a triplet, 
r = An + 2, yields n quartets and a doublet and the case r = An + 3 groups into n quartets and a triplet. Thus within 
the framework of a new kind of short-lengthed r-string (|l7|), we can provide a possible explanation of the non-Bethe 
string structures that have been observed for long but finite chains (l^Jl^,|l^]. For the Bethe string on the other 
hand, each root generally has two partners as its nearest neighbors, each of which in turn has its own two partners. 
Therefore the Bethe string can not close into any shorter substrings. 

For the alternative r-string ( |l7| ) the total momentum and the excitation energy may be given by p/ a6 ' ) = 
Z~2j=iPj, with cot if = ^> ancl 
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E (aa) = = y- 2 9j (18) 

" A ° ^1 o + - 9 i 

where the factor gj = g\ — -{j — 1), with gi = f for odd and 51 = | — ^- for even values of r. 

We note here that though the eigenfunction (Q) for solution (|l^) does not vanish exponentially in all directions and 

therefore does not give fully bound states, it should nevertheless be normalizable, since this solution can be considered 

only for long but finite chains with periodic boundary conditions (see [^5j). We believe that, the energy ( [l8| ) of our 

string is always lower than the energy of r-free magnons with rapidity Ao: rE\ == . 2 '' 1 . Though we could not prove 

0+3 

this conjecture in general, we checked it for a number of cases. Few of such sample cases are r = 4, 5, 7, 9, 10 with 
A = 1.0, whichgives E ( r as) = 1.94,1.92,2.9,3.9,4.8 and rE 1 = 3.2,4.0,5.6,7.2,8.0, correspondingly. Any new string 
solutions like (|17|) or non-Bethe excitations observed earlier, if they exist, should possibly arise from the interactions 
of different Bethe strings at the thermodynamic limit, since the later form a complete set. 

Precise formulation of the valid Bethe string ( |l2| ) with (|l3|), (|lj) and presentation of a new string ( |l7j) are the main 
results of this paper. Solution ( |l7j ) can be generalized to strings with nonuniform spacing, since at the thermody- 
namic limit it splits functionally into smaller groups with not much relevance on their relative distances. Non-Bethe 
macroscopic string with nonuniform spacing has been considered recently [pof . The results obtained here for the 
simplest case of XXX-spm chain can be extended also to other integrable models allowing Bethe string solutions, 
like XXZ-spm chain, ID Hubbard model, supersymmetric t — J model etc. and field models like NLS, SG, Liouville 
model. Therefore, application of the TBA by using the alternative string proposed here to the models investigated 
earlier with the use of Bethe strings Q would be an interesting problem. 
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